3.309 \(\int \frac{x^2 \tanh ^{-1}(a x)^2}{(1-a^2 x^2)^3} \, dx\)

Optimal. Leaf size=163 \[ -\frac{x}{64 a^2 \left (1-a^2 x^2\right )}+\frac{x}{32 a^2 \left (1-a^2 x^2\right )^2}-\frac{x \tanh ^{-1}(a x)^2}{8 a^2 \left (1-a^2 x^2\right )}+\frac{x \tanh ^{-1}(a x)^2}{4 a^2 \left (1-a^2 x^2\right )^2}+\frac{\tanh ^{-1}(a x)}{8 a^3 \left (1-a^2 x^2\right )}-\frac{\tanh ^{-1}(a x)}{8 a^3 \left (1-a^2 x^2\right )^2}-\frac{\tanh ^{-1}(a x)^3}{24 a^3}-\frac{\tanh ^{-1}(a x)}{64 a^3} \]

[Out]

x/(32*a^2*(1 - a^2*x^2)^2) - x/(64*a^2*(1 - a^2*x^2)) - ArcTanh[a*x]/(64*a^3) - ArcTanh[a*x]/(8*a^3*(1 - a^2*x
^2)^2) + ArcTanh[a*x]/(8*a^3*(1 - a^2*x^2)) + (x*ArcTanh[a*x]^2)/(4*a^2*(1 - a^2*x^2)^2) - (x*ArcTanh[a*x]^2)/
(8*a^2*(1 - a^2*x^2)) - ArcTanh[a*x]^3/(24*a^3)

________________________________________________________________________________________

Rubi [A]  time = 0.250392, antiderivative size = 163, normalized size of antiderivative = 1., number of steps used = 13, number of rules used = 6, integrand size = 22, \(\frac{\text{number of rules}}{\text{integrand size}}\) = 0.273, Rules used = {6028, 5956, 5994, 199, 206, 5964} \[ -\frac{x}{64 a^2 \left (1-a^2 x^2\right )}+\frac{x}{32 a^2 \left (1-a^2 x^2\right )^2}-\frac{x \tanh ^{-1}(a x)^2}{8 a^2 \left (1-a^2 x^2\right )}+\frac{x \tanh ^{-1}(a x)^2}{4 a^2 \left (1-a^2 x^2\right )^2}+\frac{\tanh ^{-1}(a x)}{8 a^3 \left (1-a^2 x^2\right )}-\frac{\tanh ^{-1}(a x)}{8 a^3 \left (1-a^2 x^2\right )^2}-\frac{\tanh ^{-1}(a x)^3}{24 a^3}-\frac{\tanh ^{-1}(a x)}{64 a^3} \]

Antiderivative was successfully verified.

[In]

Int[(x^2*ArcTanh[a*x]^2)/(1 - a^2*x^2)^3,x]

[Out]

x/(32*a^2*(1 - a^2*x^2)^2) - x/(64*a^2*(1 - a^2*x^2)) - ArcTanh[a*x]/(64*a^3) - ArcTanh[a*x]/(8*a^3*(1 - a^2*x
^2)^2) + ArcTanh[a*x]/(8*a^3*(1 - a^2*x^2)) + (x*ArcTanh[a*x]^2)/(4*a^2*(1 - a^2*x^2)^2) - (x*ArcTanh[a*x]^2)/
(8*a^2*(1 - a^2*x^2)) - ArcTanh[a*x]^3/(24*a^3)

Rule 6028

Int[((a_.) + ArcTanh[(c_.)*(x_)]*(b_.))^(p_.)*(x_)^(m_)*((d_) + (e_.)*(x_)^2)^(q_), x_Symbol] :> Dist[1/e, Int
[x^(m - 2)*(d + e*x^2)^(q + 1)*(a + b*ArcTanh[c*x])^p, x], x] - Dist[d/e, Int[x^(m - 2)*(d + e*x^2)^q*(a + b*A
rcTanh[c*x])^p, x], x] /; FreeQ[{a, b, c, d, e}, x] && EqQ[c^2*d + e, 0] && IntegersQ[p, 2*q] && LtQ[q, -1] &&
IGtQ[m, 1] && NeQ[p, -1]

Rule 5956

Int[((a_.) + ArcTanh[(c_.)*(x_)]*(b_.))^(p_.)/((d_) + (e_.)*(x_)^2)^2, x_Symbol] :> Simp[(x*(a + b*ArcTanh[c*x
])^p)/(2*d*(d + e*x^2)), x] + (-Dist[(b*c*p)/2, Int[(x*(a + b*ArcTanh[c*x])^(p - 1))/(d + e*x^2)^2, x], x] + S
imp[(a + b*ArcTanh[c*x])^(p + 1)/(2*b*c*d^2*(p + 1)), x]) /; FreeQ[{a, b, c, d, e}, x] && EqQ[c^2*d + e, 0] &&
GtQ[p, 0]

Rule 5994

Int[((a_.) + ArcTanh[(c_.)*(x_)]*(b_.))^(p_.)*(x_)*((d_) + (e_.)*(x_)^2)^(q_.), x_Symbol] :> Simp[((d + e*x^2)
^(q + 1)*(a + b*ArcTanh[c*x])^p)/(2*e*(q + 1)), x] + Dist[(b*p)/(2*c*(q + 1)), Int[(d + e*x^2)^q*(a + b*ArcTan
h[c*x])^(p - 1), x], x] /; FreeQ[{a, b, c, d, e, q}, x] && EqQ[c^2*d + e, 0] && GtQ[p, 0] && NeQ[q, -1]

Rule 199

Int[((a_) + (b_.)*(x_)^(n_))^(p_), x_Symbol] :> -Simp[(x*(a + b*x^n)^(p + 1))/(a*n*(p + 1)), x] + Dist[(n*(p +
1) + 1)/(a*n*(p + 1)), Int[(a + b*x^n)^(p + 1), x], x] /; FreeQ[{a, b}, x] && IGtQ[n, 0] && LtQ[p, -1] && (In
tegerQ[2*p] || (n == 2 && IntegerQ[4*p]) || (n == 2 && IntegerQ[3*p]) || Denominator[p + 1/n] < Denominator[p]
)

Rule 206

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(1*ArcTanh[(Rt[-b, 2]*x)/Rt[a, 2]])/(Rt[a, 2]*Rt[-b, 2]), x]
/; FreeQ[{a, b}, x] && NegQ[a/b] && (GtQ[a, 0] || LtQ[b, 0])

Rule 5964

Int[((a_.) + ArcTanh[(c_.)*(x_)]*(b_.))^(p_)*((d_) + (e_.)*(x_)^2)^(q_), x_Symbol] :> -Simp[(b*p*(d + e*x^2)^(
q + 1)*(a + b*ArcTanh[c*x])^(p - 1))/(4*c*d*(q + 1)^2), x] + (Dist[(2*q + 3)/(2*d*(q + 1)), Int[(d + e*x^2)^(q
+ 1)*(a + b*ArcTanh[c*x])^p, x], x] + Dist[(b^2*p*(p - 1))/(4*(q + 1)^2), Int[(d + e*x^2)^q*(a + b*ArcTanh[c*
x])^(p - 2), x], x] - Simp[(x*(d + e*x^2)^(q + 1)*(a + b*ArcTanh[c*x])^p)/(2*d*(q + 1)), x]) /; FreeQ[{a, b, c
, d, e}, x] && EqQ[c^2*d + e, 0] && LtQ[q, -1] && GtQ[p, 1] && NeQ[q, -3/2]

Rubi steps

\begin{align*} \int \frac{x^2 \tanh ^{-1}(a x)^2}{\left (1-a^2 x^2\right )^3} \, dx &=\frac{\int \frac{\tanh ^{-1}(a x)^2}{\left (1-a^2 x^2\right )^3} \, dx}{a^2}-\frac{\int \frac{\tanh ^{-1}(a x)^2}{\left (1-a^2 x^2\right )^2} \, dx}{a^2}\\ &=-\frac{\tanh ^{-1}(a x)}{8 a^3 \left (1-a^2 x^2\right )^2}+\frac{x \tanh ^{-1}(a x)^2}{4 a^2 \left (1-a^2 x^2\right )^2}-\frac{x \tanh ^{-1}(a x)^2}{2 a^2 \left (1-a^2 x^2\right )}-\frac{\tanh ^{-1}(a x)^3}{6 a^3}+\frac{\int \frac{1}{\left (1-a^2 x^2\right )^3} \, dx}{8 a^2}+\frac{3 \int \frac{\tanh ^{-1}(a x)^2}{\left (1-a^2 x^2\right )^2} \, dx}{4 a^2}+\frac{\int \frac{x \tanh ^{-1}(a x)}{\left (1-a^2 x^2\right )^2} \, dx}{a}\\ &=\frac{x}{32 a^2 \left (1-a^2 x^2\right )^2}-\frac{\tanh ^{-1}(a x)}{8 a^3 \left (1-a^2 x^2\right )^2}+\frac{\tanh ^{-1}(a x)}{2 a^3 \left (1-a^2 x^2\right )}+\frac{x \tanh ^{-1}(a x)^2}{4 a^2 \left (1-a^2 x^2\right )^2}-\frac{x \tanh ^{-1}(a x)^2}{8 a^2 \left (1-a^2 x^2\right )}-\frac{\tanh ^{-1}(a x)^3}{24 a^3}+\frac{3 \int \frac{1}{\left (1-a^2 x^2\right )^2} \, dx}{32 a^2}-\frac{\int \frac{1}{\left (1-a^2 x^2\right )^2} \, dx}{2 a^2}-\frac{3 \int \frac{x \tanh ^{-1}(a x)}{\left (1-a^2 x^2\right )^2} \, dx}{4 a}\\ &=\frac{x}{32 a^2 \left (1-a^2 x^2\right )^2}-\frac{13 x}{64 a^2 \left (1-a^2 x^2\right )}-\frac{\tanh ^{-1}(a x)}{8 a^3 \left (1-a^2 x^2\right )^2}+\frac{\tanh ^{-1}(a x)}{8 a^3 \left (1-a^2 x^2\right )}+\frac{x \tanh ^{-1}(a x)^2}{4 a^2 \left (1-a^2 x^2\right )^2}-\frac{x \tanh ^{-1}(a x)^2}{8 a^2 \left (1-a^2 x^2\right )}-\frac{\tanh ^{-1}(a x)^3}{24 a^3}+\frac{3 \int \frac{1}{1-a^2 x^2} \, dx}{64 a^2}-\frac{\int \frac{1}{1-a^2 x^2} \, dx}{4 a^2}+\frac{3 \int \frac{1}{\left (1-a^2 x^2\right )^2} \, dx}{8 a^2}\\ &=\frac{x}{32 a^2 \left (1-a^2 x^2\right )^2}-\frac{x}{64 a^2 \left (1-a^2 x^2\right )}-\frac{13 \tanh ^{-1}(a x)}{64 a^3}-\frac{\tanh ^{-1}(a x)}{8 a^3 \left (1-a^2 x^2\right )^2}+\frac{\tanh ^{-1}(a x)}{8 a^3 \left (1-a^2 x^2\right )}+\frac{x \tanh ^{-1}(a x)^2}{4 a^2 \left (1-a^2 x^2\right )^2}-\frac{x \tanh ^{-1}(a x)^2}{8 a^2 \left (1-a^2 x^2\right )}-\frac{\tanh ^{-1}(a x)^3}{24 a^3}+\frac{3 \int \frac{1}{1-a^2 x^2} \, dx}{16 a^2}\\ &=\frac{x}{32 a^2 \left (1-a^2 x^2\right )^2}-\frac{x}{64 a^2 \left (1-a^2 x^2\right )}-\frac{\tanh ^{-1}(a x)}{64 a^3}-\frac{\tanh ^{-1}(a x)}{8 a^3 \left (1-a^2 x^2\right )^2}+\frac{\tanh ^{-1}(a x)}{8 a^3 \left (1-a^2 x^2\right )}+\frac{x \tanh ^{-1}(a x)^2}{4 a^2 \left (1-a^2 x^2\right )^2}-\frac{x \tanh ^{-1}(a x)^2}{8 a^2 \left (1-a^2 x^2\right )}-\frac{\tanh ^{-1}(a x)^3}{24 a^3}\\ \end{align*}

Mathematica [A]  time = 0.116072, size = 121, normalized size = 0.74 \[ \frac{6 a x \left (a^2 x^2+1\right )+3 \left (a^2 x^2-1\right )^2 \log (1-a x)-3 \left (a^2 x^2-1\right )^2 \log (a x+1)-16 \left (a^2 x^2-1\right )^2 \tanh ^{-1}(a x)^3+48 \left (a^3 x^3+a x\right ) \tanh ^{-1}(a x)^2-48 a^2 x^2 \tanh ^{-1}(a x)}{384 a^3 \left (a^2 x^2-1\right )^2} \]

Antiderivative was successfully verified.

[In]

Integrate[(x^2*ArcTanh[a*x]^2)/(1 - a^2*x^2)^3,x]

[Out]

(6*a*x*(1 + a^2*x^2) - 48*a^2*x^2*ArcTanh[a*x] + 48*(a*x + a^3*x^3)*ArcTanh[a*x]^2 - 16*(-1 + a^2*x^2)^2*ArcTa
nh[a*x]^3 + 3*(-1 + a^2*x^2)^2*Log[1 - a*x] - 3*(-1 + a^2*x^2)^2*Log[1 + a*x])/(384*a^3*(-1 + a^2*x^2)^2)

________________________________________________________________________________________

Maple [C]  time = 0.441, size = 2571, normalized size = 15.8 \begin{align*} \text{result too large to display} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

int(x^2*arctanh(a*x)^2/(-a^2*x^2+1)^3,x)

[Out]

1/64/(a*x-1)^2/(a*x+1)^2*x^3+1/16/a^3*arctanh(a*x)^2/(a*x-1)^2-1/16/a^3*arctanh(a*x)^2/(a*x+1)^2+1/32*I*a/(a*x
-1)^2/(a*x+1)^2*arctanh(a*x)^2*Pi*csgn(I*(a*x+1)/(-a^2*x^2+1)^(1/2))^2*csgn(I*(a*x+1)^2/(a^2*x^2-1))*x^4+1/16*
I*a/(a*x-1)^2/(a*x+1)^2*arctanh(a*x)^2*Pi*csgn(I*(a*x+1)/(-a^2*x^2+1)^(1/2))*csgn(I*(a*x+1)^2/(a^2*x^2-1))^2*x
^4+1/32*I*a/(a*x-1)^2/(a*x+1)^2*arctanh(a*x)^2*Pi*csgn(I*(a*x+1)^2/(a^2*x^2-1))^3*x^4+1/32*I*a/(a*x-1)^2/(a*x+
1)^2*arctanh(a*x)^2*Pi*csgn(I*(a*x+1)^2/(a^2*x^2-1)/((a*x+1)^2/(-a^2*x^2+1)+1))^3*x^4-1/16*I*a/(a*x-1)^2/(a*x+
1)^2*arctanh(a*x)^2*Pi*csgn(I/((a*x+1)^2/(-a^2*x^2+1)+1))^3*x^4+1/16*I*a/(a*x-1)^2/(a*x+1)^2*arctanh(a*x)^2*Pi
*csgn(I/((a*x+1)^2/(-a^2*x^2+1)+1))^2*x^4-1/16*I/a/(a*x-1)^2/(a*x+1)^2*arctanh(a*x)^2*Pi*csgn(I*(a*x+1)^2/(a^2
*x^2-1))^3*x^2-1/16*I/a/(a*x-1)^2/(a*x+1)^2*arctanh(a*x)^2*Pi*csgn(I*(a*x+1)^2/(a^2*x^2-1)/((a*x+1)^2/(-a^2*x^
2+1)+1))^3*x^2+1/8*I/a/(a*x-1)^2/(a*x+1)^2*arctanh(a*x)^2*Pi*csgn(I/((a*x+1)^2/(-a^2*x^2+1)+1))^3*x^2-1/8*I/a/
(a*x-1)^2/(a*x+1)^2*arctanh(a*x)^2*Pi*csgn(I/((a*x+1)^2/(-a^2*x^2+1)+1))^2*x^2+1/32*I/a^3/(a*x-1)^2/(a*x+1)^2*
Pi*arctanh(a*x)^2*csgn(I*(a*x+1)/(-a^2*x^2+1)^(1/2))^2*csgn(I*(a*x+1)^2/(a^2*x^2-1))+1/16*I/a^3/(a*x-1)^2/(a*x
+1)^2*Pi*arctanh(a*x)^2*csgn(I*(a*x+1)/(-a^2*x^2+1)^(1/2))*csgn(I*(a*x+1)^2/(a^2*x^2-1))^2-1/32*I/a^3/(a*x-1)^
2/(a*x+1)^2*Pi*arctanh(a*x)^2*csgn(I*(a*x+1)^2/(a^2*x^2-1))*csgn(I*(a*x+1)^2/(a^2*x^2-1)/((a*x+1)^2/(-a^2*x^2+
1)+1))^2+1/32*I/a^3/(a*x-1)^2/(a*x+1)^2*Pi*arctanh(a*x)^2*csgn(I/((a*x+1)^2/(-a^2*x^2+1)+1))*csgn(I*(a*x+1)^2/
(a^2*x^2-1)/((a*x+1)^2/(-a^2*x^2+1)+1))^2+1/16/a^3*arctanh(a*x)^2/(a*x-1)+1/16/a^3*arctanh(a*x)^2*ln(a*x-1)+1/
16/a^3*arctanh(a*x)^2/(a*x+1)-1/16/a^3*arctanh(a*x)^2*ln(a*x+1)+1/8/a^3*arctanh(a*x)^2*ln((a*x+1)/(-a^2*x^2+1)
^(1/2))-1/16*I/a^3/(a*x-1)^2/(a*x+1)^2*Pi*arctanh(a*x)^2-1/32*I*a/(a*x-1)^2/(a*x+1)^2*arctanh(a*x)^2*Pi*csgn(I
*(a*x+1)^2/(a^2*x^2-1))*csgn(I*(a*x+1)^2/(a^2*x^2-1)/((a*x+1)^2/(-a^2*x^2+1)+1))^2*x^4+1/32*I*a/(a*x-1)^2/(a*x
+1)^2*arctanh(a*x)^2*Pi*csgn(I/((a*x+1)^2/(-a^2*x^2+1)+1))*csgn(I*(a*x+1)^2/(a^2*x^2-1)/((a*x+1)^2/(-a^2*x^2+1
)+1))^2*x^4-1/16*I/a/(a*x-1)^2/(a*x+1)^2*arctanh(a*x)^2*Pi*csgn(I*(a*x+1)/(-a^2*x^2+1)^(1/2))^2*csgn(I*(a*x+1)
^2/(a^2*x^2-1))*x^2-1/8*I/a/(a*x-1)^2/(a*x+1)^2*arctanh(a*x)^2*Pi*csgn(I*(a*x+1)/(-a^2*x^2+1)^(1/2))*csgn(I*(a
*x+1)^2/(a^2*x^2-1))^2*x^2+1/16*I/a/(a*x-1)^2/(a*x+1)^2*arctanh(a*x)^2*Pi*csgn(I*(a*x+1)^2/(a^2*x^2-1))*csgn(I
*(a*x+1)^2/(a^2*x^2-1)/((a*x+1)^2/(-a^2*x^2+1)+1))^2*x^2-1/16*I/a/(a*x-1)^2/(a*x+1)^2*arctanh(a*x)^2*Pi*csgn(I
/((a*x+1)^2/(-a^2*x^2+1)+1))*csgn(I*(a*x+1)^2/(a^2*x^2-1)/((a*x+1)^2/(-a^2*x^2+1)+1))^2*x^2-1/32*I/a^3/(a*x-1)
^2/(a*x+1)^2*Pi*arctanh(a*x)^2*csgn(I/((a*x+1)^2/(-a^2*x^2+1)+1))*csgn(I*(a*x+1)^2/(a^2*x^2-1))*csgn(I*(a*x+1)
^2/(a^2*x^2-1)/((a*x+1)^2/(-a^2*x^2+1)+1))+1/64/a^2/(a*x-1)^2/(a*x+1)^2*x-1/24/a^3/(a*x-1)^2/(a*x+1)^2*arctanh
(a*x)^3-1/64/a^3/(a*x-1)^2/(a*x+1)^2*arctanh(a*x)-1/32*I*a/(a*x-1)^2/(a*x+1)^2*arctanh(a*x)^2*Pi*csgn(I/((a*x+
1)^2/(-a^2*x^2+1)+1))*csgn(I*(a*x+1)^2/(a^2*x^2-1))*csgn(I*(a*x+1)^2/(a^2*x^2-1)/((a*x+1)^2/(-a^2*x^2+1)+1))*x
^4+1/16*I/a/(a*x-1)^2/(a*x+1)^2*arctanh(a*x)^2*Pi*csgn(I/((a*x+1)^2/(-a^2*x^2+1)+1))*csgn(I*(a*x+1)^2/(a^2*x^2
-1))*csgn(I*(a*x+1)^2/(a^2*x^2-1)/((a*x+1)^2/(-a^2*x^2+1)+1))*x^2-1/24*a/(a*x-1)^2/(a*x+1)^2*arctanh(a*x)^3*x^
4-1/64*a/(a*x-1)^2/(a*x+1)^2*arctanh(a*x)*x^4+1/12/a/(a*x-1)^2/(a*x+1)^2*arctanh(a*x)^3*x^2-3/32/a/(a*x-1)^2/(
a*x+1)^2*arctanh(a*x)*x^2+1/32*I/a^3/(a*x-1)^2/(a*x+1)^2*Pi*arctanh(a*x)^2*csgn(I*(a*x+1)^2/(a^2*x^2-1))^3+1/3
2*I/a^3/(a*x-1)^2/(a*x+1)^2*Pi*arctanh(a*x)^2*csgn(I*(a*x+1)^2/(a^2*x^2-1)/((a*x+1)^2/(-a^2*x^2+1)+1))^3-1/16*
I/a^3/(a*x-1)^2/(a*x+1)^2*Pi*arctanh(a*x)^2*csgn(I/((a*x+1)^2/(-a^2*x^2+1)+1))^3+1/16*I/a^3/(a*x-1)^2/(a*x+1)^
2*Pi*arctanh(a*x)^2*csgn(I/((a*x+1)^2/(-a^2*x^2+1)+1))^2-1/16*I*a/(a*x-1)^2/(a*x+1)^2*arctanh(a*x)^2*Pi*x^4+1/
8*I/a/(a*x-1)^2/(a*x+1)^2*arctanh(a*x)^2*Pi*x^2

________________________________________________________________________________________

Maxima [B]  time = 1.00358, size = 524, normalized size = 3.21 \begin{align*} \frac{1}{16} \,{\left (\frac{2 \,{\left (a^{2} x^{3} + x\right )}}{a^{6} x^{4} - 2 \, a^{4} x^{2} + a^{2}} - \frac{\log \left (a x + 1\right )}{a^{3}} + \frac{\log \left (a x - 1\right )}{a^{3}}\right )} \operatorname{artanh}\left (a x\right )^{2} + \frac{{\left (6 \, a^{3} x^{3} - 2 \,{\left (a^{4} x^{4} - 2 \, a^{2} x^{2} + 1\right )} \log \left (a x + 1\right )^{3} + 6 \,{\left (a^{4} x^{4} - 2 \, a^{2} x^{2} + 1\right )} \log \left (a x + 1\right )^{2} \log \left (a x - 1\right ) + 2 \,{\left (a^{4} x^{4} - 2 \, a^{2} x^{2} + 1\right )} \log \left (a x - 1\right )^{3} + 6 \, a x - 3 \,{\left (a^{4} x^{4} - 2 \, a^{2} x^{2} + 2 \,{\left (a^{4} x^{4} - 2 \, a^{2} x^{2} + 1\right )} \log \left (a x - 1\right )^{2} + 1\right )} \log \left (a x + 1\right ) + 3 \,{\left (a^{4} x^{4} - 2 \, a^{2} x^{2} + 1\right )} \log \left (a x - 1\right )\right )} a^{2}}{384 \,{\left (a^{9} x^{4} - 2 \, a^{7} x^{2} + a^{5}\right )}} - \frac{{\left (4 \, a^{2} x^{2} -{\left (a^{4} x^{4} - 2 \, a^{2} x^{2} + 1\right )} \log \left (a x + 1\right )^{2} + 2 \,{\left (a^{4} x^{4} - 2 \, a^{2} x^{2} + 1\right )} \log \left (a x + 1\right ) \log \left (a x - 1\right ) -{\left (a^{4} x^{4} - 2 \, a^{2} x^{2} + 1\right )} \log \left (a x - 1\right )^{2}\right )} a \operatorname{artanh}\left (a x\right )}{32 \,{\left (a^{8} x^{4} - 2 \, a^{6} x^{2} + a^{4}\right )}} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(x^2*arctanh(a*x)^2/(-a^2*x^2+1)^3,x, algorithm="maxima")

[Out]

1/16*(2*(a^2*x^3 + x)/(a^6*x^4 - 2*a^4*x^2 + a^2) - log(a*x + 1)/a^3 + log(a*x - 1)/a^3)*arctanh(a*x)^2 + 1/38
4*(6*a^3*x^3 - 2*(a^4*x^4 - 2*a^2*x^2 + 1)*log(a*x + 1)^3 + 6*(a^4*x^4 - 2*a^2*x^2 + 1)*log(a*x + 1)^2*log(a*x
- 1) + 2*(a^4*x^4 - 2*a^2*x^2 + 1)*log(a*x - 1)^3 + 6*a*x - 3*(a^4*x^4 - 2*a^2*x^2 + 2*(a^4*x^4 - 2*a^2*x^2 +
1)*log(a*x - 1)^2 + 1)*log(a*x + 1) + 3*(a^4*x^4 - 2*a^2*x^2 + 1)*log(a*x - 1))*a^2/(a^9*x^4 - 2*a^7*x^2 + a^
5) - 1/32*(4*a^2*x^2 - (a^4*x^4 - 2*a^2*x^2 + 1)*log(a*x + 1)^2 + 2*(a^4*x^4 - 2*a^2*x^2 + 1)*log(a*x + 1)*log
(a*x - 1) - (a^4*x^4 - 2*a^2*x^2 + 1)*log(a*x - 1)^2)*a*arctanh(a*x)/(a^8*x^4 - 2*a^6*x^2 + a^4)

________________________________________________________________________________________

Fricas [A]  time = 2.05374, size = 294, normalized size = 1.8 \begin{align*} \frac{6 \, a^{3} x^{3} - 2 \,{\left (a^{4} x^{4} - 2 \, a^{2} x^{2} + 1\right )} \log \left (-\frac{a x + 1}{a x - 1}\right )^{3} + 12 \,{\left (a^{3} x^{3} + a x\right )} \log \left (-\frac{a x + 1}{a x - 1}\right )^{2} + 6 \, a x - 3 \,{\left (a^{4} x^{4} + 6 \, a^{2} x^{2} + 1\right )} \log \left (-\frac{a x + 1}{a x - 1}\right )}{384 \,{\left (a^{7} x^{4} - 2 \, a^{5} x^{2} + a^{3}\right )}} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(x^2*arctanh(a*x)^2/(-a^2*x^2+1)^3,x, algorithm="fricas")

[Out]

1/384*(6*a^3*x^3 - 2*(a^4*x^4 - 2*a^2*x^2 + 1)*log(-(a*x + 1)/(a*x - 1))^3 + 12*(a^3*x^3 + a*x)*log(-(a*x + 1)
/(a*x - 1))^2 + 6*a*x - 3*(a^4*x^4 + 6*a^2*x^2 + 1)*log(-(a*x + 1)/(a*x - 1)))/(a^7*x^4 - 2*a^5*x^2 + a^3)

________________________________________________________________________________________

Sympy [F]  time = 0., size = 0, normalized size = 0. \begin{align*} - \int \frac{x^{2} \operatorname{atanh}^{2}{\left (a x \right )}}{a^{6} x^{6} - 3 a^{4} x^{4} + 3 a^{2} x^{2} - 1}\, dx \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(x**2*atanh(a*x)**2/(-a**2*x**2+1)**3,x)

[Out]

-Integral(x**2*atanh(a*x)**2/(a**6*x**6 - 3*a**4*x**4 + 3*a**2*x**2 - 1), x)

________________________________________________________________________________________

Giac [F]  time = 0., size = 0, normalized size = 0. \begin{align*} \int -\frac{x^{2} \operatorname{artanh}\left (a x\right )^{2}}{{\left (a^{2} x^{2} - 1\right )}^{3}}\,{d x} \end{align*}

Verification of antiderivative is not currently implemented for this CAS.

[In]

integrate(x^2*arctanh(a*x)^2/(-a^2*x^2+1)^3,x, algorithm="giac")

[Out]

integrate(-x^2*arctanh(a*x)^2/(a^2*x^2 - 1)^3, x)